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ABSTRACT 

Numerical simulation and analysis of the solu- 
tion are presented for a ;aminar reacting flowfield of 
air and hydrogen in the case of an external combus- 
tion employed to reduce base drag in hypersonic ve- 
hicles operating at transonic speeds. The flowfield 
consists of a transonic air stream at a Mach num- 
ber of 1.26, and a sonic transverse hydrogen injection 
along a row of 26 orifices. Self-sustained combustion 
is computed over an expansion ramp downstream of 
the injection and a flaneholder, using the recently 
developed RPLUS code. Measured data is available 
only for surface pressure distributions, and is used 
for validation of the coae in practical 3D reacting 
flowfields. Pressure comparison shows generally good 
agreements and the mam effects of combustion are 
also qualitatively consistent with experiment. 


INTRODUCTION 

Results obtained from a three-dimensional com- 
putation of hydrogen-air combustion are presented 
for flowfields simulating an external burning exper- 
iment tl], which was performed to study external 
combustion as a mechanism for drag reduction in 
aerospace vehicles when operating in the transonic 
regime. The paper presents an analysis of the numer- 
ical solutions computed to validate the recently de- 
veloped RPLUS code. The feasibility of the proposed 
concept for drag reduction is not of primary concern 
in this paper, but qualitative conclusions however can 
be drawn from the results for certain prominent ef- 
fects of combustion. 

One of the important components in hypersonic 
vehicles is the nozzle which must be designed to oper- 
ate in a wide range of flow regimes. At the on-design 
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hypersonic conditions, considerably large area ratios 
are required to accommodate for the extremely high 
pressure ratios needed to obtain rapid expansion to 
hypersonic Mach numbers [2,3]. Due to weight con- 
sideration of such area ratio requirement, the aft-end 
of the aircraft is utilized as the expansion surface 
for the upper nozzle wall. The other lower surface 
must also be truncated, producing a shear layer be- 
tween an internal plume exhaust and an external am- 
bient flow. In the transonic regime which occurs at 
low speed off-design conditions, the nozzle becomes 
severely overexpanded [3], leading to negative thrust 
or nozzle drag. One possible means to reduce the in- 
ternal overexpansion is to use variable nozzle geome- 
try by deflecting the truncated surface. Although the 
nozzle drag then can be reduced, a base drag, how- 
ever, is created as a consequence of the ambient flow 
expansion on the external side of the truncated sur- 
face. External combustion is, therefore, proposed as 
a mechanism to remove this base drag and to affect 
favorably the exhaust plume in order to achieve fully 
expanded flowfields. In ref. [1], an analysis of the 
potential impact of the external burning concept has 
been presented, and an experiment has also been con- 
ducted to study combustion if it is attainable under 
the transonic conditions. 

The experimental configuration consists simply 
of a small segment of a flat plate connected to a 12- 
degree expansion ramp, see fig. (1). Over these three- 
dimensional surfaces, an air stream flows at a super- 
sonic Mach number of 1.26 measured at the leading 
edge. A sonic transverse injection of hydrogen fuel 
is located along a row of 26 orifices placed in front 
of the expansion corner. Downstream of the injec- 
tion on the flat plate, a rectangular block to be used 
as flameholder is installed with the trailing edge po- 
sitioned at the expansion corner. The flameholder 
has a small opening at the symmetry plane necessary 
for thermal expansion. The described geometry can 
be viewed as a combined representation of an exter- 
nal undersurface of the truncated nozzle wall and a 
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shear layer occurring between the exhaust and ambi- 
ent flows. 

The overall flowfield contains a leading edge 
shock wave and a downstream Prandtl-Meyer expan- 
sion fan in the absence of combustion. When there 
is combustion, the expansion fan is weaken or elimi- 
nated, which is the desired effect making the external 
air stream to flow without turning at the corner. 

At a Mach number of 1.26, the leading edge 
shock normally detaches because of the boundary 
layer. More importantly, due to the experimental 
setup, the underside of the model is a compression 
surface also inclined at a 12-degree angle. The sur- 
face then produces a high pressure region below the 
model, causing further shock detachment from the 
leading edge. Consequently, a considerably strong 
upwash at the ieading edge is observed in the exper- 
iment. The incoming stream therefore is not parallel 
to the flat plate but flows locally at an angle of attack 
near the leading edge. Simulating this effect is costly 
for the full moael since the computational domain 
becomes rather large. A uniform supersonic profile 
at Mach 1.26 is then assumed here, for simplicity, as 
inflow condition at the leading edge. 

The RPLUS code [4] is employed to calculate 
the results presented in this paper. Development and 
validation of the code are still progressing. Currently, 
turbulence models have not been implemented, and 
therefore only laminar flows can be calculated. How- 
ever, significant physical features are reasonably cap- 
tured in the present computations although turbu- 
lence can be incorporated in future work to improve 
the results quantitatively. The code has a capabil- 
ity of multi-block grids for treating more complicated 
surface boundaries, such as in regions of the flame- 
holder or the tnermal expansion slot in the present 
geometry where single grids are difficult to generate. 
To obtain the solution, the code utilizes an LU solver 
and takes advantage of 2D vectorization on an oblique 
grid plane of the computational domain. Although 
the vectorization is coded only for the implicit side of 
the solver, the overall procedure can lead to an effi- 
cient alternative for calculating, particularly, reacting 
flows. 

Some of the recent works on complex 3D re- 
acting flows using the RPLUS code have been doc- 
umented for supersonic combustors [5] and for a hy- 
drogen fueled ramjet engine [6], where the flows occur 
generally in tne supersonic regime. This paper illus- 
trates another validation case through an application 
of the code to practical reacting flowfields in the tran- 


sonic regime. 


NUMERICAL METHOD 

The RPLUS code solves for the solutions of the 
full compressible Navier-Stokes and species eolations 
in a curvilinear coordinate system. A finite /olume 
scheme is formulated in the discretization he gov- 
erning equations. The method is secon rder ac- 

curate on a uniform grid, equivalent to u ^g central 
differences in finite difference formulations. Jameson- 
type artificial dissipations must also be added to 
avoid numerical oscillation and instability. For the 
solutions, an iterative LU algorithm is employed, re- 
quiring only scalar inversions of the flow equations 
but block diagonal inversions of the species equations 
because of the chemical source terms. All variables 
are dimensional. 

Finite-rate chemistry of hydrogen and air is 
based on a combustion model containing 18 reaction 
steps involving 9 species with nitrogen assumed in- 
ert. The gas mixture is taken to be thermally perfect. 
Specific heat, thermal conductivity and viscosity are 
curve-fitted by fourth degree polynomials of temper- 
ature for each species. Mixture specific heat is evalu- 
ated by species mass concentration weighting, while 
mixture thermal conductivity and viscosity are ob- 
tained from Wilke’s mixing rule. Binary diffusion co- 
efficient is computed based on Chapman- Enskog the- 
ory using Lennard-Jones potential functions. More 
detailed description of the RPLUS code can be found 
in ref. [4j. 


GEOMETRY AND GRID 

Figure (1) schematically illustrates the exper- 
imental geometry and the corresponding system of 
coordinates. The model has a total lengtn of 36.07 
cm and a width of 15.24 cm. The flat plate section 
is 7.62 cm in length measured from the leading edge, 
and the following ramp is inclined at an angle of 12 
degrees. The fiameholder with a trailing edge at the 
corner juncture has a narrow length of 0.635 cm and 
a height of 0.318 cm, and it is disjointed in the middle 
by the small thermal expansion slot having 0.318 cm 
in width. Hydrogen fuel is sonically injected along a 
row of 26 circular orifices, which is located at 1.27 cm 
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in front of the expansion corner. These orifices are 
uniformly distributed across the model width, each 
having a diameter of 0.064 cm. The experimental 
pressure probes are positioned along the symmetry 
plane. 

The corresponding numerical simulation has the 
inflow and outflow boundaries which coincide with 
the leading edge of the flat plate and the trailing 
edge of the expansion ramp respectively. Shock de- 
tachment is therefore not allowed since the flow in 
front of the leading edge is not included in the com- 
putational domain. Similarly, the wake downstream 
of the expansion ramp is also not simulated. 

Only one half of the model needs to be com- 
puted in the z direction, and one of the boundaries 
in this direction is the symmetry plane. The compu- 
tational domain taken in the present work includes 
the region exterior to the model side edge. The sec- 
ond boundary is then a farfield uniform flow. In 
the y direction, the upper boundary is located at 
a distance sufficiently far from the surface so that 
an undisturbed freestream can be prescribed. The 
lower boundary coincides with the bottom edge of 
the model. 

A grid representative of the desenoed compu- 
tational domain is depicted in fig. ( 21, illustrat- 
ing surface and boundary grid distributions. Only 
parts of the inflow, symmetry and lower boundaries 
are shown, for clarity. Large grid density across the 
model width is required to resolve the injection ori- 
fices. Similarly in the x direction, grid concentration 
occurs at the injection and around the flameholder. A 
close-up view of the region around the expansion cor- 
ner is presented in fig. (3), showing the grids for the 
rectangular- block flameholder, the thermal expansion 
opening between one end of the flameholder and the 
symmetry plane, and a series of the injection orifices 
indicated by the very small shaded surfaces. The ex- 
perimental circular orifice is numerically modeled as 
a square which has the same surface area to obtain 
the same injected mass flow rate. Each orifice is ap- 
proximated by two grid cells in the x direction and 
one cell in the z direction. For gradual distribution, 
several grid points are also needed in the space be- 
tween the injectors. A large grid size is therefore re- 
quired to resolve all 13 orifices for one half of the full 
model, and consequently the computation time be- 
comes considerable for reacting flows. Although grid 
refinement for the injectors is possible, no attempt is 
made here in this regard since demand for computa- 
tion time with such further refinement is exceedingly 


high. 

Two types of multi-block grids have been used 
for the computations. One is the internal blocking 
which is inherent in the RPLUS code and is formu- 
lated for complex geometries where the computation 
domain is divided into several blocks defined by the 
surface contours. This type of blocking is used to 
simulate the flameholder and the thermal expansion 
slot. The other is the external blocking by separat- 
ing the zones of combustion, which takes place near 
the ramp surface, and nonreacting flow, which occurs 
above the combustion zone and extends into the in- 
viscid region. The partitioning is simple since grid 
coordinates between the two zones are matched iden- 
tically, although mismatched grid can also be con- 
structed and has been used in ref. [5j. Presently, the 
use of matched grids allows straightforward transfer- 
ring of solutions at the interface without the need for 
interpolation. Interaction between the two zones is 
iteratively performed external to RPLUS. The main 
purpose of external blocking is to reduce computer 
time particularly in the nonreacting region since there 
is only one species, air, to be calculated and conver- 
gence can be reached much faster than in the com- 
bustion zone. External blocking is also utilized here 
in order to place the upper boundary at a distance 
far above the ramp surface with coarse mesh distri- 
butions. 

The total grid has a dimension of 48x50x96, 
with the largest dimension in the z direction hav- 
ing 84 grid points distributed across the model width. 
External blocking is used only in the y direction, with 
10 grid points in the nonreacting zone near the up- 
per boundary, and 30 points in the combustion zone 
above the ramp surface. The region exterior to the 
side edge has 10 grid points distributed from the lower 
boundary to the ramp surface. The smallest spacing 
is about 0.0001 cm at the no-slip surfaces, sufficient 
to resolve the laminar viscous layer. 


RESULTS 

Two cases of reacting and nonreacting flowfields 
have been computed for the same flow' conditions on 
an identical grid. The incoming air is provided at a 
Mach number of 1.26, a temperature of 228 K and 
a pressure of 0.54 atm. whereas hydrogen fuel is in- 
jected at a Mach of 1, a temperature of 228 K and a 
pressure of 23.8 atm. Because of the low temperature 
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used in both streams, spontaneous chemical reaction 
of air and hydrogen cannot occur, but the injection 
leads to a mixing, nonreacting fiowfield which does 
not deviate significantly from the case of pure air 
flows except in the region near the injection. For self- 
sustained combustion to take place, some form of heat 
source must be introduced initially to induce chemical 
reaction. Experimentally, spark ignition was used but 
possible only at low subsonic incoming Mach num- 
bers. The air stream was then gradually adjusted to 
the supersonic condition of Mach 1.26. For the react- 
ing results presented here, combustion is numerically 
initiated using higher fuel temperature and then re- 
ducing it to the desired temperature of 228 K after 
combustion has already established in the fiowfield. 
An ignition fuel temperature value of 1400 K has been 
used successfully to achieve the initial induced reac- 
tion and the subsequent self-sustain combustion. Re- 
sults of the mixing and reacting cases are presented 
and compared in the discussion below. 

Figure (4) illustrates a typical streamwise flow 
pattern in terms of Mach contours in a z-constant 
plane between the side edge and the symmetry plane 
for the nonreacting case. Only a portion of the flow- 
field near the surface is shown and discussed herein. 
The basic structure contains a very weak leading-edge 
shock completely diffused by numerical dissipations, 
and a well-defined expansion fan at the corner. Ef- 
fects of hydrogen injection are negligible in the in- 
viscid region above the injection, where it consists 
only of pure air. A smali region of separation is 
also visible behind the flameholder at this particu- 
lar spanwise (z) location. Fig. (5) depicts a three- 
dimensional view for Mach contours at various cross- 
sectional planes. Only subsonic Mach numbers are 
plotted so that comparisons can be made latter with 
the combustion case. It can be observed that the 
flow is highly three-dimensional with the existence 
of streamwise vortices indicated by two large regions 
of subsonic flows. The vortex along the symmetry 
plane is produced because of the thermal expansion 
slot, while the other, smaller vortex is present due to 
the geometry and the pressure difference between the 
ramp surface and the ambient flow. 

When there is combustion, the fiowfield has a 
significantly different structure. Mach contours for 
the reacting case are shown on the same z-plane in 
fig. (6), and at various x-planes in fig. (7) plotted, 
again, only for subsonic Macn numbers. Importantly, 
both of the vortices are diminished and the flow shows 
two-dimensionality within :ne model. The inviscid 


region above the ramp is markedly altered, containing 
a weaken expansion fan having a smaller turning an- 
gle, and a subsonic zone which increases substantially 
in size and becomes more uniform across the model 
width. However, upstream influence of combustion 
appears to be minimal in front of the flameholder. 

The effect of combustion, which is to prevent 
flow turning at the expansion corner, can be seen by 
comparing fig. (8) and (9) plotted for the particle 
traces in the same z-piane used for Mach numbers 
above. First, it can be seen clearly that the stream- 
line for the noncombustion case is somewhat parallel 
to the ramp surface because of flow turning. With 
combustion, the streamline remains nearly horizon- 
tal behind the flameholder and it is then shifted up- 
wards in a region near the outflow where combus- 
tion is more extensive. Second, combustion enlarges 
the separated region, as seen in fig. (9). The en- 
largement of the separated flow also extends in the 
spanwise direction which can be seen from the next 
figures (10,11), In the nonreacting case, fig. (10), 
the separation is centered near the symmetry plane, 
whereas in the reacting case, fig. (11), the bubble is 
centered near the side edge, showing separated flows 
occurring nearly throughout the model width. The 
streamwise vortex along the side edge is also appar- 
ent for the nonreacting flow in fig. (10), evident by 
the collapsing of the panicle traces along the model 
edge. 

Contours for hydrogen mass fraction without 
combustion are presented in fig. (12). The contour 
variation shows a fuel distribution which is uniform 
near the flameholder and then concentrated down- 
stream along the symmetry plane due to the three- 
dimensional relief effect of the high pressure ambient 
flow. Typical hydrogen mass-fraction distributions in 
the streamwise direction for the nonreacting and re- 
acting cases are shown in figs. (13,14) for the same 
z-plane midway between the symmetry plane and the 
side edge. Identical contour levels are used in both 
figures. A sparse concentration of contours far down- 
stream of the flameholder in fig. (14) indicates hy- 
drogen consumption by combustion. Top views of the 
hydrogen distribution are shown in figs. (15,16) for 
both cases, showing the relief effect without combus- 
tion and a more uniform spanwise variation in the 
combustion case. 

Contours of temperature and mass fraction of 
water vapor are illustrated in figs. (17,18). Both fig- 
ures depict a similar pattern where large variations 
exist across a thick layer above the ramp surface, 
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where high rates of reaction occur. In general, The 
flow reaches a peak temperature in a region near the 
outflow and slightly above the ramp surface. Typi- 
cal top view of temperature contours is presented in 
fig. (19) at a vertical distance near the ramp surface. 
The dense distribution along the symmetry plane cor- 
responds to a narrow region originated from the ther- 
mal expansion slot, in which the combustion is very 
weak largely due to high concentration of hydrogen 
with little mixing. Combustion is more extensive in 
the interior region having a peak temperature occur- 
ring near the side edge at the outflow. Fig. (20) 
demonstrates the contours for mass fraction of water 
vapor in several cross-sectional planes. Combustion 
intensity increases for the flow in regions near the 
wall where highest concentration of water vapor is 
produced. 

Figure (21) depicts distributions of mass flow 
rates for hydrogen ana some of the species whose pro- 
ductions are relatively high. For hydrogen, compar- 
ison is made between reacting and nonreacting. As 
can be seen, hydrogen mass flow rate remains con- 
stant without combustion, while gradual reduction 
of hydrogen is evident with combustion. Approxi- 
mately, forty percent of hydrogen is consumed in the 
case of a laminar flow assumed in the present com- 
putation. Water vapor is predominantly present in 
the combustion products but, for other species, the 
productions are very small or negligible. The cor- 
responding distributions for mass fraction of chemi- 
cal species and temperature are shown in fig. (22). 
The variations represent the maximum values taken 
at each axial computational plane and then plot- 
ted along the ramp length. Mass fraction for wa- 
ter vapor is also prominent in magnitude as com- 
pared to those of other species, with a largest value 
located near the flameholder and a decreasing vari- 
ation downstream. Axially maximum temperature 
increases toward the outflow, showing sudden com- 
bustion behind the flameholder and continuing reac- 
tion accumulated as the flow progresses. The peak 
temperature of about 2050 K is in agreement with 
experiment. 

Temperature distributions are also plotted in 
figs. (23.24) for several z-planes from the symmetry 
plane at 2=0, to near the model side edge at z=0.075. 
The plotted values are also locally maximum, taken 
along vertical grid lines at each axial plane. At z=0 as 
in fig. (23), it can be clearly seen that the combustion 
takes place more gradually than abruptly, consisting 
of two subsequent Temperature rises with strong re- 


action occurring latter in the second stage initiated 
from the middle section of the ramp. For the interior 
region up to z=0.055 in fig. (24), the distributions 
are similar to one another and combustion is rather 
sudden behind the flameholder. The presence of the 
flameholder is, therefore, to keep the flame front to at- 
tach at the expansion corner rather than downstream 
as exhibited at the symmetry plane. Regions asso- 
ciated with high temperature are located near the 
side edge as indicated by the curve at z=0.064 in fig. 
(24). It is noted that this curve is similar to that 
in fig. (22). Next to this region, the distributions 
at z=0.071 and 0.075 represent the variations at a 
location corresponding to the spanwise boundary of 
the combustion zone along the side edge where the 
reaction starts to diminish. 

Comparison with experiment is made in figs. 
(25,26) for surface pressure distributions along the 
symmetry plane for both nonreacting and reacting 
cases. The computed results show qualitatively cor- 
rect trends of the pressure variations according to 
measurement. The discrepancy near the leading edge 
is attributed primarily to the existence of an upwash 
at the leading edge where experimental data shows 
a negative pressure coefficient. Apparently, the ex- 
perimental flow has a local angle of attack at the 
leading edge, which is critically high to cause a very 
large pressure drop. The use of a uniform super- 
sonic profile as inflow condition to simplify the sim- 
ulation also contributes to the discrepancy since ef- 
fect of shock detachment is not taken into account. 
There is also disagreement behind the expansion cor- 
ner, and it is reasonably attributed to the laminar 
regime which can lead to severe regions of separa- 
tion associated with adverse pressure gradients. In 
general, the agreement between computation and ex- 
periment is good, considering the flow complexity and 
the simplifying approximations employed in the simu- 
lation. Most importantly, the computed results pre- 
dict fairly well the effect of combustion, that is to 
increase pressure on the ramp surface, as can be seen 
by comparing figs. (25) and (26) particularly near 
the corner. Note that the peak pressure corresponds 
to the location of hydrogen injection. Measured pres- 
sure in the combustion case remains constant behind 
the flameholder. whereas computed pressure shows 
oscillation due to the large region of laminar sepa- 
ration. Surface pressure tends to approach ambient 
pressure near the outflow. For nonreacting, this re- 
covery is caused by the three-dimensional relief effect 
previously mentioned. Combustion in the experiment 
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also exhibits considerable upstream influence up to 
the leading edge, suggesting that, on the flat plate, 
there is a thick subsonic region adjacent to the viscous 
layer or even a reversed flow. This subsonic flow per- 
mitting upstream influence is created by the detached 
bow shock and further affected by the upwash, but 
is not present in the computed flowfield since a very 
thin subsonic layer is produced because of a uniform 
supersonic inflow. 

The presented results for the nonreacting case 
have been converged to a residual reduction of nearly 
three orders of magnitude. Further reduction is possi- 
ble but the convergence rate becomes slow. Each iter- 
ation takes 14 seconds and the total time is about 100 
hours on a YMP computer. For the reacting case, the 
convergence has been quite difficult and extremely 
slow. Each iteration requires 45 seconds and the to- 
tal time is approximately 400 hours for a reduction 
of only two orders of magnitude. Under- relaxation 
of the temperature and pressure fields is used. The 
lack of robustness and consequently the slow conver- 
gence rate are likely caused by the chemical source 
terms, even though partially implicit treatments of 
these terms in the RPLUS code have reduced some 
of the inherent stiffness. Another possibility worth 
mentioning is the outflow conditions. In the reacting 
case, the subsonic region at the outflow is substantial, 
and the boundary condition here is not determined 
unless the wake behind the ramp must also be simu- 
lated. Extrapolation as used in the present work can 
adversely affect the convergence behavior. 
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SUMMARY 


Numerical results obtained from the RPLUS 
code have been presented for a case of external com- 
bustion using air and hydrogen. Both nonreacting 
and reacting flows in the laminar regime are simu- 
lated and surface pressure is compared to experimen- 
tal data for validation. The agreement between mea- 
sured and computed pressure is qualitatively good. 
Most importantly, the effect of combustion, that is to 
increase the ramp pressure and consequently to pre- 
vent flow turning at the expansion juncture, is pre- 
dicted consistently with experiment. 
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fig. 4 Mach contours, z-plane, nonreacting. 
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fig. 6 Mach contours, z-plane, reacting. 
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fig, 25 Pressure distribution, symmetry plane, nonreacting. 
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